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Abstract 

We review some aspects of the theory of spherical Bessel functions and Struve functions by means 
of an operational procedure essentially of umbral nature, capable of providing the straightforward 
evaluation of their definite integrals and of successive derivatives. The method we propose allows 
indeed the formal reduction of these family of functions to elementary ones of Gaussian type. We 
study the problem in general terms and present a formalism capable of providing a unifying point 
of view including Anger and Weber functions too. The link to the multi-index Bessel functions is 
also briefly discussed. 
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I. INTRODUCTION 



The spherical Bessel and Struve functions [L2J are widely exploited in applications, e. 
g. the diffraction and scattering of radiation 3j. Albeit their properties are rather well 
mown, we apply here a new point of view, developed in a recent series of papers (see Ref. 
4(), which frames them within the operational context. This method is fairly attractive and 
yields a significant simplification of the underlying computations. 

In the first part of the paper we deal with the theory of spherical Bessel functions, while 
the second will be devoted to the study of the properties of the Struve functions. The 
underlying thread will be discussed in the concluding section. 

The spherical Bessel functions j n (x) are linked to their cylindrical counterparts by jl, 2, 5| 



3n(x) 




J n+ i (x) 
2x +2 



in e Z). 



The relevant differential equation can easily be derived as follows. We start by introducing 
the operators f| 

E± = — =F d x , 

x 

where N is an operator whose action on a generic indexed function is given by iy G M) 



NF v {x) = vF v {x). 



By using the relation [l| 

J' v {x) = =f 
we obtain the following identity 



Jv±i{x) Ju{x) 

X 



E± JJx) = J v ±\{x) 



1.21 



and, therefore 



E+E-JJx) = JJx) 



In the case v = n+1/2, taking into account of definition (11.11) . this equation can be rewritten 
in differential terms as follows 



U ~\ Q \ ( n+ \ 



X 



x 



d x 



<2x 



7T 



■3n{X) 



<2x 



7l 



■3n[X) 
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and, after a few manipulations, we end up with 

{x 2 d 2 x + 2xd x + [x 2 -n(n + l)]} j n (x) = 0. 
According to the definition (II. ip . the function j n (x) can be expressed by the series 

JnW " V2x ^Jfc!r(n + Jb + |) V2. 
which can be formally written as 



with the operator c defined by identity 

Although we have already exploited this method for other families of Bessel functions Q], 
we get a glimpse of how the formal re-handling of the spherical Bessel functions given in Eq. 
(II. 3p may be useful, by considering the evaluation of the integral 



jo(x)dx. 

According to Eq. (11.31) . by treating c as an ordinary constant we can reduce our problem to 
the evaluation of a Gaussian integral, obtaining 

oo 

e -e(*/2) 2 dx 



V ,(0) = ^^^(0) = 7 rcV(0) = 7 r. 
In the forthcoming sections we will further prove the usefulness of this procedure. 

II. DERIVATIVES AND INTEGRALS OF SPHERICAL BESSEL FUNCTIONS 

The operational definition (II. 3p of j n (x) allows a very straightforward derivation of the 
relevant generating function. We get 

oo 



t n V 7T C f C 1 

^ n\ 3n{x) = ~F exp ~i (x2 " 2xt) r ^ (0) 

1 =-L = jo ( Vx 2 -2xt) . 2.1 

2 Vx 2 -2xt J \ ) 



n=0 

' 7T 
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The formalism of the shift operators Q] can now be usefully exploited to get further prop- 
erties. By recalling indeed that 

exp | A (±d z \ | f(z) = f (v^+2A) , 

Eq. ( 12. ip can be rewritten as 



n=0 ' n=0 



( -d x ) j (x) 



ml \x 



from which we obtain the well-known property 



i=xt 



jn(x) = (~x) n ( -d x J j (x). 

The described approach also allows us to easily obtain a closed form for the successive 
derivative of the function jo(x). By taking into account the following identity 9| 



d n x e ax = H n (2ax,a) e aQ 



where H n (y, z) are the two- variable Hermite polynomials defined by 

fc=0 v ' 

from Eq. (II. 3p we get 



d n x Jo(x) 



fc!(n-2A:)! Jn " H J ' 



c(a;/2) 2 



v(o) 



As a final example of application of the method, let us now consider the following integral 

/oo 
j n (x)dx. 
-oo 

By using the generating function method we can write 



n=0 



71 C 



exp ^—-(x — 2 x t) \ dx 



<p(0) 



oo ^ / ^ \ 2 k 



TV 



fc=0 



(fc! 



(2.2) 
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from which, by equating the coefficients of the same powers in t, we get 
bo n = V-^r ^r- = -V T n + - , 6, n+ , = 0. 



From Eq. (12 .2p we also get 

b(t)=1il (t) 

where I (t) is the modified Bessel function of order and, thus, as a consequence of Eq. 
(12. ip . the following identity results 

7T _ / b 



jo (^Jax 2 + bx^j dx 



/n 



la \2 ■s/a / 

that can also be viewed as an integral representation of the function Jo- 



in. THE STRUVE FUNCTIONS 

The Struve functions H a (x) are defined by the series (see Eq. (57:6:1) in Ref. 

(-l) k 



n a ( x ) = 



rj- \ 2 m+a+1 

2. 



fe= -r(m+|) r(m + a + |) 

The derivation of the relevant differential equation can be achieved through the shift oper- 
ators defined in Eq. ( 11. 2(1 . By using the differentiation formula (see Eq. (57:10:1) in Ref. 

) 

(x/2)° 



d x HJx) 



1 

2 



Ha-Ax) - H 



a+l 



X + 



^T(a+|) 

and taking into account the recursion formula for the Struve functions 



2 a 

H a+1 (x) + H a _i(x) = — H Q (x) + 

x 



(x/2Y 



r(a + |)' 



it is easy to show that 

£■+ H a (x) = H a+1 (x) 



(x/2Y 



(« + |)' 



£_ H Q (x) = Ha_i(l) 



i.e. 



E. E_ H n (x) = HJx) 



[x/2) 



a-l 



In differential terms the last equation can be rewritten as follows 



a 



0., 



x 



a-l 



d x ) H Q (x) = H Q (x) - 



(x/2) 



-X 

2> 
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that can be reduced to the following non-homogeneous Bessel equation 

4 (x/2) a+1 



B a HJx) 



(x d x ) + (x 



2 a 2 ) 



v^r(a + i)' 

The properties of Struve functions can be studied by means of a slight modification of 
the procedure we have followed so far. To widen the perspective we consider the problem 
starting from the study of the Humbert-Bessel functions , which are two- index Bessel-like 



functions defined by the series 



k=0 



-X) 



k\T(k + (i + l)T(k + v+l) 



(3.1) 



and whose connection with the Struve functions is realized by the following identity 



H.(x) = (I)"" jf 



e~ s Ji n , i 

2J L 2' a +2 



x 



ds. 



According to our technique, the Humbert-Bessel functions can formally written as 

1 



J a,v\X) — Ci Co e 



IJ, ~V — Cl C2 x 



Pi (0)^(0), c?<pi(0) 



1,2), (3.2) 



r(l + a) v 

an expression that makes a very simple task to derive, for example, the identities reported 
below 



J^ v (x )dx 



r (/, + !) r(„ + i)' 



io r(/x — a + 1) T(z/ — a + 1) 

The operatorial expression (13.21) allows one also to easily prove that (see Eq. (2.7.2.1) in 



Ref. 



ul) 



m n (x) dx 



1 + I) c~ {a+1)/2 ct m j o e~° s- 1 -*' 2 ds 



cot a 



-2 < Re(a) < 0), 



4j, where are also 



which can be further handled to get the canonical form reported in Ref. 
discussed the conditions of validity. 

A more general definition of Struve-like functions can be obtained through the definition 
of the following auxiliary function 

.,2- 



/ e s 1 Jafi ys — ) ds 



X 
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i.e., by using Eq. (I3.2p 



sJ -i e - s( i +£l c 2X yA) ds 



¥>1 (0)^(0) 



r(7) 



X 



X F 2 ( 7; 1 + a, 1 + # -— } . 



r(i + a)r(i + /3) 

with ((a p ); (6 9 ); 2) the generalized hypergeometric function. From the definition (13. ip it 
is also easy to show that 



00 

V w m v n J m , n (x) =exp{w + t;- — ) 



m,n=— 00 



and, therefore, the generating function of A aii g )7 (x) with respect to the indices a and /3 is 
given by 

1 /X\ 2 



^ />oo r 

u m v n A m>n>1 {x) = e u+v / exp -1 

r(7) 



m,n=— 00 



7 — s 



ii« V2 



ds 



u+v 



7 ■ 



uv \2J 

The multi-index Bessel of Humbert type functions are not widespreadly known, although 
their use could be very advantageous in applications. As an example, we consider the case 
of the product of two cylindrical Bessel functions, which can be written as (see Eq. 8.442.1 



from 



5|) 



J„{x) J v {x) = 



■l) k T(fi + u + 2k + l) 



k=0 



k\r(fi + k + l)T(u + k + l)r(fi + u + k + l) \2J 



i.e. 



with 



-S l±+U T / 



* 2X 1 ds 



(3.3) 



E 

fc=0 



{-zf 



k\ T{fi + k + 1) I> + k + 1) r(p + k + 1) ' 



In the context of our operatorial method, this series can be written as 

•W*) = K^K^ cx ^ iZ <pi(o) ^2(0)^3(0) 

that, inserted in Eq. fl 3 . 3 j) . makes an easy task to prove, for example, the following identity 



x\ -v-v 



Jn(x) J„(x)dx = \pK 



T(n + u) 



r u + - r [V + -] r u + z/ + - 
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IV. CONCLUDING REMARKS 



Multi-index Bessel-like functions can also be used to express the Anger J u (x) and Weber 



Ej,(x) functions, defined as 



10| 




cos 



sin 



v x 
~2~ 

V X 

IT 



sin 



cos 



v x 
v x 




with 



Si{v,x) = 



l) k (x/2) 



2k 



IT Ik 



v 
2 



S 2 {u,x) = Y 



-l) k (x/2) 2k+1 



^=0 r [ k + 



r k 



V 

2 



These functions can be expressed in the operatorial form (Si^ (i, k 
symbol) 



1,2) is the Kronecker 



S k {v, x) 



A5 k ,2+v)/2 A5 k>2 -u)/2 



(x/2) 



8k;. 



^(0)^(0) 



1 + C! c 2 (x/2) 2 

that can be used to easily prove, for example, the following identities 

S 2 {viX 



J Si(u,x)dx = cos J , y 



dx = — sin ( — — 



.r 



7T 



In this paper we have shown that the use of concepts of umbral nature from the opera- 
tional calculus, combined with the properties of special functions and polynomials, can be 
successfully applied to the the theory of Bessel-like functions. The method we have sug- 
gested provides some advantages in "practical" computations and is flexible enough to open 
many new perspectives, which cannot all be explored in the space of a single paper. The 
main drawback is the lack of mathematical rigor and, therefore, the necessity of checking 
the obtained results resorting, for example, to numerical methods. 

Some of our results (those relevant to the integrals) can be justified on the basis of the 
Ramanujan master theorem 12Hl4|. originally suggested by XIX century operationalists like 
Glaisher p|. However, we believe that the strategy one should follow, when dealing with 
the present method, is that of using it as a procedure to first "guess" some specific formulae 
and then use them as sound conjecture to be proven with other conventional means. 
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